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Abstract. Yang-Mills theory in four dimensions is studied by using the Coulomb gauge. 
The Coulomb gauge Hamiltonian involves integration of matrix elements of an operator 



> 

O 

Q ■ V built from the Laplacian and from a first-order differential operator. The operator 

O i V is studied from the point of view of spectral theory of pseudo-differential operators 



on compact Riemannian manifolds, both when self-adjointness holds and when it is not 
fulfilled. In both cases, well-defined matrix elements of V are eval 
towards the more difficult problems of quantized Yang-Mills theory. 



^' fulfilled. In both cases, well-defined matrix elements of V are evaluated as a first step 

> 

X 



1. Introduction 

Although much progress has been made in the theoretical understanding of gauge theories 
of fundamental interactions/"^ a number of outstanding unsolved problems remain. In 
particular, the mass gap problem has been stressed in the recent literature,^ which amounts 
to proving that, for any compact simple gauge group, quantum Yang-Mills theory on R"^ 
has an Hamiltonian operator with no spectrum in the interval (0, 6) for some 6 > 0. One 
should then show that, starting from the classical action functional 



-^ I TriF^^F^d'x, (1) 



the corresponding Hamiltonian operator has energy spectrum bounded from below, with 
strictly positive lower bound. The solution of such a problem is very important not only 
to achieve internal consistency of the mathematical formalism, but also (if not mainly) 
for physical reasons: the existence of a mass gap would account for the nuclear force 
being strong but short-ranged. Moreover, the recent theory of glueballs^ relies entirely 
on action functionals like (1), without fermionic fields. In plain terms, it is possible for 
two or more gluons to combine into a strongly bound, neutral-coloured particle of pure 
glue. This (hypothetical) object is called a glueball.^ Moreover, a gluon can bind with a 
meson to form a hybrid. For example, a red quark and an anti-blue antiquark can bind 
with a blue/antired gluon. The lightest glueball allowed by quantum chromo dynamics can 
be described by a circular tube of glue and has vanishing angular momentum.^ This has 
spherical symmetry while glueballs of other, elongated shapes have no n- vanishing angular 
momenta and larger masses. Hybrids can announce their presence by yielding at least three 
s-wave mesons. They are in fact predicted to decay into one s-wave meson and another 
short-lived meson with internal angular momentum. The latter then breaks up into two s- 
wave mesons. The first (though uncertain) experimental evidence in favour of hybrids was 
obtained in 1994, when experimenters at Protvino found an object called 7r(1800), emerging 
from collisions of pions with protons. This particle has the quantum properties and decay 
pattern expected for a hybrid. Recent research deals with special hybrids called exotics, 
which have combinations of internal angular momentum, parity and charge conjugation 



quantum numbers that are forbidden for mesons. For example, the simplest exotic has 
J = 1,P = — 1 and C = 1. Inconclusive evidence also exists in favour of glueballs called 
/o(1500) and /j(1710), which would belong to the class of glueballs of mass in the range 
1.500 to 1.800 MeV.^ The latest experiments tell us that glueball candidates and qq mesons 
have been found to be produced with different momentum and angular dependences in the 
central region of pp collisions.^ 

All this phenomenology can be described, in principle, with the help of a Yang-Mills 
Lagrangian. We have been therefore motivated, in our research, by the mass gap prob- 
lem in four dimensions, although it will become clear that, for the time being, we only 
have a possibly new perspective in classical Yang-Mills theory. For this purpose, we have 
considered the Yang-Mills Lagrangian in the Coulomb gauge with the associated decom- 
position into electric and magnetic parts. Although the Coulomb gauge is non-covariant, 
it has the merit of leading to a quantum theory which is manifestly unitary. Moreover, 
well-established calculational recipes are available for writing down the Hamiltonian oper- 
ator of the quantum theory,^ while recent work has shown that the Coulomb gauge can be 
viewed as the singular limit of the Landau-Coulomb interpolating gauge, ^ adding therefore 
evidence in favour of such a gauge being very appropriate for the quantum theory. 

Section 2 presents a review of the Coulomb gauge Hamiltonian, with emphasis on ma- 
trix elements of an operator V which plays a key role in Sec. 3. Here, attention is focused on 
compact Euclidean spacetime, for which discrete spectral resolutions of self-adjoint elliptic 
operators exist. The matrix elements ofV among square-integrable functions are evaluated 
explicitly in such a case. Section 4 presents an assessment of the spectral approach to the 
mass gap. 

2. Coulomb Gauge Hamiltonian 

To help the general reader, we present here a brief review of the Coulomb gauge Hamilto- 
nian in classical and quantum theory, relying on Ref. 7. 



The Lagrangian L of a classical Yang-Mills field can be expressed in terms of electric 
and magnetic fields, which are defined by 

8, = -—A, - {V^ + gA,x)Ao = -—A, + -V,u, (2) 

at at g 

EjkiBi = VjAk - VkAj + glj X At. (3) 

In Eq. (2), Aq = — -, and Vi is the covariant derivative in the Coulomb gauge: 

V, = V^ + gA,x, (4) 

the cross denoting the isovector product in all our equations. Thus, the Lagrangian reads 

L^]^j{s,-8,-B,-B^d\, (5) 

with corresponding Hamiltonian 

H = ^f(s,-S, + B,-B,y^r. (6) 

One can now exploit the constraint VjAj = and decompose the electric field into a 
transverse part and a gradient term, i.e. 

S, = Sr - V,0, (7) 

where (hereafter we define the Laplacian as the operator A = — V'^Vfc, with a minus sign 
in front of second derivatives to make it bounded from below) 

£l' = -ll,+ (<5,, + A-iV,V,)(l, X 00), (8) 

which satisfies ViS'^^ = 0, and 

= lo - <?A-i (1, X V,Ao) . (9) 



Since the electric field has vanishing divergence with respect to the connection D, i.e. 
ViSi = 0, one gets also an equation for 0, i.e. 



9 



(10) 



Such an equation is solved upon inverting the operator V^Pi, which yields (summation 
over repeated indices is understood) 



giVkVur^A, X £f^ = -g{VkVk)-^A, x nf . 



:iii 



At this stage our analysis is still classical, but in the quantum theory the counterpart of 
Eq. (11) is more involved because zero- modes of the Faddeev-Popov operator play a role 
and their effect should be included. 

Upon exploiting the transverse nature of £\^ and imposing fall-off conditions at infinity 
on (h, the classical Hamiltonian reads 



// = ! 

2 



e + V 



2 / ^\2 



d^r. 



(12) 



This is further simplified by using the identity 



(V,0)(V,(/)) = V,((/)V,0) + A 0, 



(13) 



jointly with Eq. (11). Thus, on defining 



n" 



(14) 



which represents the charge carried by Aj, the classical Hamiltonian reads eventually 



2 



nr +[B, 



d^r 



+ Y I (y\{r){l,r\V\l\r')a^A{r')d^rd^r' . 



(15) 



In Eq. (15), V is the integro-differential operator defined by 

V = {V,V,)-^ /\{VPj)-\ (16) 

where the inverse of V^Pi results from (11), and the Laplacian results from (13). The 
scheme obtained from (15) is not amenable to calculation unless one finds a convenient 
way of expressing the inverse of 'Wi'Di. For this purpose, one defines 

A = V,(i'^x) = Ai X Vi, (17) 

so that 

V{D, = - A +^A, (18) 

and hence 

(V.A)-' = (- A +gK)-^ = -A-i (/ + ^A A'^ H9^/\-^? + 0(/)) 

- - A-i -g A-i A A-i -g'^ A'^ A A'^ A A'^ +0{g^), (19) 

where we have applied the formula 

{AB)-^ = B-^A-^ 

to the operators A ^ I — gAA~^ and B = —A. The insertion of the (formal) expansion 
(19) into the definition (16) yields an algorithm for V, i.e. 

oo 

P = A-i + ^(A; + l)/A-i (AA-l)^ (20) 

which is useful at small g and for the analysis of spectral asymptotics. 

In the quantum theory, one has instead to consider the un-renormalized coupling 
constant go and the Faddeev-Popov determinant^ 

7 = det(V.A). (21) 



The equal-time canonical commutation relations read (the smeared form is more rigorous 
but inessential for our purposes, which are not axiomatic) 



A, (f, ty, Ul'if', tr] = t5'^ ((5,fc + A-iy, Vfc) d^{r- r'), 



(22) 



A{f,ty,A,{r',tr = ^f(f,t)^^f(r-^,t) 



and the Hamiltonian operator in the Coulomb gauge takes the form^ 



(23) 



H 



-Ifrtr 



^ 



^J-'W-lW + ^B:. 



d^r 



+ ^ f-f-^a\r}{l,r\V\l',r)W'{r)d^rd^r. 



(24) 



Note that, since quantum fields are operator-valued distributions,^ the Hamiltonian as in 
the form just written is not defined as an operator in a Hilbert space. Strictly, the product 
of local operators should be regularized, and all quantum formulae should be written with 
this understanding. 

3. Structure of V 

Our ultimate goal is the investigation of the spectrum of H, with the associated role played 
by 7, n*'^, Bi,aK But this is still extremely difficult, and hence we resort to the analysis of 
the operator V (see Eq. (16)) occurring in the classical theory, here written in the form 



V = {A-gK)-^/\{A-gK)-\ 



(25) 



where one should bear in mind that A is the first-order differential operator defined in Eq. 
(17). Now two main cases can be distinguished in a framework where Minkowski space-time 
is replaced by a compact four- geometry (M, g) without boundary. This is more relevant 
for Euclidean field theory, which is nevertheless an important branch of quantum field 
theory. ^'^'^^ By doing so, one may hopefully learn lessons about operators whose analysis 
is a mandatory step before being able to solve the original problem. 



(i) Assume first that the operator V is self-adjoint. Since, under the above assumptions on 
{M,g), the Laplacian A is self-adjoint, this means we are treating the isovector product 
Ai X Vi as a self-adjoint operator (see definition (17)), so that 

V^ = (At - ^At)-i A^ (At - ^At)-i = V. (26) 

One can then exploit theorems ensuring that the eigenfunctions ui of V form a complete or- 
thonormal set. Any square-integrable function (f G L'^{M) can be then expanded according 
to (q being the Fourier coefficients q = {ui, ip)) 

oo 

ip = '^ciui. (27) 

1=1 

The resulting mean value of V reads (A^ being its eigenvalues, for which Vui = XiUi) 

oo oo 

{f,Vif) = J2 >^iCrCi{ur,ui) = ^A^IqI^, (28) 

l,r=l 1=1 

while more general matrix elements read (with bi = {ui, $)) 

oo 

($,Pv^) = 5^Az6rQ, (29) 

1=1 

for all $ and ip G LF'{M). The mean value of V is therefore positive if its spectrum is 
bounded from below, with positive lower bound. 

(ii) Even when {M,g) is compact, the operator A — gA may fail to be self-adjoint, since 
A contains the effect of first-order covariant derivatives, which may be anti-self-adjoint. If 
this were the case, we can nevertheless exploit the expansion (20) to point out that V has 
the general structure 

V = A-^+V\ (30) 

where V is a pseudo-differential operator (see Appendix) of order —3, since the lowest 
value of k in the infinite sum in Eq. (20) involves the operator A~^ A A~^. One can then 



rely upon the work of Beals/^ who has studied spectral asymptotics of elliptic operators 
with a self-adjoint positive principal part, showing that the eigenvalues lie in a parabolic 
region around the positive axis. In our case the principal part of V is the inverse of the 
Laplacian, which is both self-adjoint and positive. Such a qualitative information should 
be supplemented by well known results about spectral resolutions of partial differential 
operators of positive order on compact manifolds. ^^ In other words, lack of self-adjointness 
of V makes it now impossible to expand Lp G LF'{M) according to Eq. (27), while an 
expansion of Lp remains legitimate if a discrete spectral resolution of the Laplacian is used. 
For this purpose, let us denote by vi the eigenfunctions of the Laplacian with discrete 
eigenvalues ni, for which 

Avi = iJiVi, (31) 



and {Ci being the Fourier coefficient {vi,(p)) 



oo 



<f^J2CiVi. (32) 

1=1 



With the help of the expansion (20), which is useful at small g, the mean value of the 
operator V is then found to take the form 



oo 

|2 



CXD oo ^ 

+ E ^rCi Y^ik + 1)/- [vr, (AA-i)SO . (33) 

l,r=l k=l ^^ 

Such a formula shows how the mean value of V changes on passing from the self-adjoint 
(see Eq. (28)) to the non-self-adjoint case. The formula (33) cannot be further simplified, 
because the operator A does not commute with the inverse of the Laplacian. The first few 
powers of AA~^ in the sum over all positive values of the integer k yield, for example, 

AA-^vi = —Avi, 



and such formulae are useful if one needs to stop the analysis of (</?, Pc/?) at terms of order 
0{g'^). If the formula (33) is used, positivity of the mean value of V amounts to proving 
that the following inequality holds: 

Y^-\Ci\^ > -Y,C:CiY.g'^-{vrAAA-^fvi). (34) 

Z=l ^^ l,r=l k=l ^^ 

Similarly, for all $ and tp G L?{M), the formula (29) for the matrix element ($,P(^) is 
now replaced by (with our notation, Bi = (vi, $)) 



oo -. 

1=1 ^^ 

oo oo ^ 

+ Y, KCi J2{k + l)9'-{vr, {AA-')'^vi). (35) 

l,r=l k=l ^"^ 

The effort of studying both self-adjoint and non-self-adjoint case for V is especially 
valuable if one looks at sub-regions of M with smooth boundary, since then the particular 
boundary conditions chosen might or might not ensure self-adjointness of V. 

4. Status of the Spectral Approach 

In our paper we start from a physical motivation to move gradually towards the framework 
of pseudo-differential operators and their relevance for the Hamiltonian of classical Yang- 
Mills theory but on a background manifold with positive-definite metric. We have proposed 
a spectral-oriented perspective on the classical foundations of an outstanding open problem 
in modern quantum field theory, and we have evaluated the well-defined matrix elements 
of the fundamental operator V in the self-adjoint and non-self-adjoint case. This step is 
important, since otherwise no hope exists of being able to evaluate the matrix elements 
occurring in the quantum theory. In particular, the use of discrete spectral resolutions of 
the Laplacian for studying matrix elements of V in the non-self-adjoint case is a simple 
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but non-trivial technical point that we have advocated. At least three major problems are 
now in sight. 

(i) What happens if the compact Riemannian four-manifold M is replaced by a Lorentzian 
four-manifold, e.g. Minkowski space-time. Such a "decompactification limit" is crucial to 
recover the analysis of the original problem we are interested in. 

(ii) What happens when the classical Hamiltonian (15) is replaced by the quantum Hamil- 
tonian (24), with \l,r) being the eigenfunctionals of the position operator. 

(iii) How to make manifest the role played by the gauge group in the quantum theory (in 
particular, the group structure of the operators A and V). 

We hope, however, that the spectral approach will at least show clearly the limits of 
what can be done to solve the mass gap problem,^' ^^~^^ and it is our aim to devote our 
efforts to understand whether this is really the case. 

Appendix 

To help the reader who is not familiar with spectral theory, we here recall some basic 
concepts. A linear partial differential operator P of order d can be written in the form 



\oi\<d 

where |a| = J2T=i ^^^ and 

with act a C°° function on R"^ for all a. The associated symbol is, by definition, 

p{x,i)= J2 ac.{x)C, {A3) 

\a\<d 

i.e. it is obtained by replacing the differential operator D^ by the monomial i^". The pair 
(x, ^) may be viewed as defining a point of the cotangent bundle of R"^, and the action of 

11 



P on the elements of the Schwarz space S of smooth complex-valued functions on R"^ of 
rapid decrease is given by^"^ 

Pf{x) ^ J e'^^-y>^p{x, i)f{y)dydi, [AA) 

where the dy = dyi...dym and d$, = d^i...d^rn orders of integration cannot be interchanged, 
since the integral is not absolutely convergent. 

Pseudo-differential operators are instead a more general class of operators whose sym- 
bol need not be a polynomial (and whose order is not necessarily positive) but has suitable 
regularity properties. More precisely, let S'^ be the set of all symbols p{x,^) such that p 
is smooth in the pair of variables (x,^) with compact x support, and for all (a, /3) there 
exist constants C^^p for which 

ii):i)fp(x,e)i<a,/3(i + iei)'-i^|, (^5) 

for some real (not necessarily positive) value of d. The associated pseudo-differential 
operator, defined on the Schwarz space and taking values in the set of smooth functions 
on R"^ with compact support, 

is defined in a way formally analogous to the previous integral formula for Pf{x). 

Ellipticity can be defined by means of a majorization obeyed by the inverse of the 
modulus of the symbol. In other words, let U be an open subset with compact closure in 
R"^, and consider an open subset Ui whose closure Ui is properly included in U : Ui C U. 
If p is a symbol of order d on U, it is said to be elliptic on Ui if there exists an open set 
U2 which contains Ui and positive constants Fi so that 

\p{x,0\-'<Fi{l + m)-' {A6) 

for |,^| > Fq and x G U2, where 



lei = \ g^'ix)Ub 



\ fc=i 
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The corresponding operator P is then said to be eUiptic. 
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